Exchange integrals and magnetization distribution in BaCu2X2C>7 (X=Ge,Si) 
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Estimating the intra- and inter-chain exchange constants in BaCu2X2C>7 (X=Ge,Si) by means 
of density-functional calculations within the local spin-density approximation (LSDA) we find the 
Ge-compound to be a more ideal realization of a one-dimensional spin chain with Dzyaloshinskii- 
Moriya interaction than its Si-counterpart. Both compounds have a comparable magnitude of 
inter-chain couplings in the range of 5 ... 10 K, but the nearest neighbor intra-chain exchange of 
the Ge-compound is nearly twice as large as for the Si one. Using the LSDA+!7 method we predict 
the detailed magnetization density distribution and especially remarkable magnetic moments at the 
oxygen sites. 

PACS numbers: 71.15.Ap, 71.15.Mb, 71.15.Rf, 71.20.Be, 75.10.Pq 



The study of low-dimensional magnetism is of broad in- 
terest, especially due to the crucial role of quantum fluc- 
tuations in those systems. To investigate quantum spin 
chains experimentally one needs low-dimensional model 
compounds. A representative example is BaCu2Ge2C>7ji 
or the isomorphous compound BaCu2Si2C>7^ knowing 
to be the most ideal non-organic spin chain compound 
exhibiting the Dzyaloshinskii-Moriya (DM)^ interac- 
tion. But why are some compounds more ideal than 
others? What are the remaining differences to an ideal 
one-dimensional compound? To answer these ques- 
tions we present here a microscopic calculation of the 
electronic structure of BaCu2X 2 C>7 (X=Ge,Si) within 
density-functional theory. We investigate the chemical 
trends in comparing the Ge- with the Si-compound, es- 
pecially with respect to the intra-chain and inter-chain 
tight-binding and exchange parameters. Furthermore, we 
calculate the magnetization density in detail. 

BaCu2Gc207 and BaCu2Si207 were first synthesized 
and their crystallographic structures were described in 
1993^ Both are composed of Cu04 non- ideal plaquettes 
(the Cu atom and the 4 O atoms are not perfectly copla- 
nar) which are arranged into corner-sharing chains along 
the c axis (see Fig. 1) but with a Cu-O-Cu bonding an- 
gle between 90 and 180 degrees. These bonding angles 
are <j) = 135° (124°) for the Gc(Si)-compound, and the 
zigzag chains allow for a DM interaction (due to the low 
crystallographic symmetry, space group Pnma). Corre- 
spondingly, both compounds were intensively studied as 
examples for the influence of the DM interaction in quan- 
tum spin chains. A weak fcrromagnctism was observed in 
the Ge-compound^ which allows the experimental deter- 
mination of the spin canting angle due to the DM inter- 
action to be 1.9°. On the other hand, BaCu2Si207 has 
attracted interest in connection with the finding of two 
consecutive spin reorientation transitionsA&ifi Also, two 
gaps in the excitation spectrum of the Si-compound had 
been found by neutron scattering^ and in antiferromag- 
netic resonance measurements ~ These two gaps arise 
due to the DM interaction and the symmetric anisotropy 



term with a partial compensation between them^ii^ Re- 
cently, it was the Ge-compound which allowed to test a 
prediction of Oshikawa and Affleck^ in good accuracy. 
The linewidth of the electron paramagnetic resonance 
signal diverges at low temperatures like AH oc 1/T 2 for 
an applied magnetic field perpendicular to the DM vector 
corresponding to the creation of a staggered field. But 
for an applied field parallel to the DM vector, when no 
staggered field can be created, the linewidth vanishes like 
AH cx Ti 




FIG. 1: (Color online) Crystal structure of BaCu2Ge207. The 
crystal structure is Pnma. The chains run along the c axis. 

Despite all these results, a microscopic calculation 
which determines in detail the differences between the 
Si- and the Ge-compound is lacking up to now. It will be 
presented here. Our microscopic calculations confirm the 
naive expectation that the difference in the bonding an- 
gles is the main reason that the Ge-compound is a better 
ID model system than the Si-one. In addition, we also 
predict small oxygen moments in these compounds which 
should be detectable by several experimental methods. 
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Our analysis of exchange couplings and model pa- 
rameters follows the lines given in Ref. ^] for Sr 2 Cu03 
and Ca2Cu03. We calculate the non- magnetic band- 
structure within the local density approximation (LDA) 
and perform a tight-binding (TB) analysis of the relevant 
band at the Fermi energy. Adding a Hubbard interaction 
term we may propose an effective one-band Hamiltonian 
which allows to estimate the different exchange terms. 
These are in reasonable agreement with total energy dif- 
ferences from the LSDA+[7 method. 

The scalar rclativistic band-structure calculations were 
performed using the full potential local orbital (FPLO) 
method^ In the FPLO method a minimum basis ap- 
proach with optimized local orbitals is employed, being at 
the same time a very accurate and an efficient numerical 
tool. For the present calculations we used the following 
basis set: Ba Ad5s5p : 6s6p, Cu 3s3p : AsAp3d, Ge 3s3p 
: AsAp ; 3d5s, Si 3s3p, and O 2s2p ; 3d. The inclusion of 
the Ba Adhshp and the Cu 3s3p semi-core states into the 
valence set was necessary to account for non-negligible 
core-core overlap, and the O 3d states were used to im- 
prove the completeness of the basis set. The site-centered 
potentials and densities were expanded in spherical har- 
monics up to l m .ax = f2. The LDA calculations were 
performed using the Perdew-Wang paramctrization^i In 
all calculations we used a fc-mesh of 4 x 4 x 4 = 64 k- 
points in the Brillouin zone which is sufficient for the 
band-structure and magnetization density calculations 
reported below. 
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FIG. 2: Band-structure of BaCu2Ge207 within a non- 
magnetic state using LDA. The momentum points in the 
Brillouin space are given in units of {it /a, it /b,n /c) by F = 
(0,0,0), X = (1,0,0), Y = (0,1,0), Z = (0,0,1), s = 
(1, 1, 0), U = (1, 0, 1), R = (1, 1, 1) and T = (0, 1, 1). V -> Z 
is the dispersion along the chain axis. 

The non-magnetic band-structure of BaCu2Ge20y is 
shown in Fig. 2. One may note a well isolated band 
at the Fermi level having 8 branches according to the 
8 Cu ions per elementary cell (Fig. 1). The main va- 
lence band starts at about -1 eV and extends up to -7 
cV. It is comprised of highly hybridized Cu 3d and O 



2p orbitals as it is common for cuprate compoundsii* 
The band at the Fermi level consists mainly of the Cu 
3d x 2_ y 2 orbital lying predominantly within the plane of 
the Cu04 plaquette (in the local plaquette coordinate 
system) and hybridizing slightly by pda bonds with the 
oxygen 2p orbitals directed towards Cu. One has to 
choose the non-magnetic solution to determine the TB 
parameters^ since the Hubbard correlation will be added 
later on by hand. Taking into account nearest neighbor 
hopping processes into all 3 cartesian directions and pro- 
viding a down-folding into the present Brillouin zone we 
may approximately describe the 8 branches of the disper- 
sion of the relevant band by 

„ ( k z c\ f k T a\ ( k,.b\ 

E k = ±2t 2 cos f -y- J ± 2i x cos ( -y- J ± 2t y cos ( -J- J . 

(1) 

The corresponding fit gives the parameters listed in the 
Table [I] The errors are estimated by the evident de- 
viations of the real band-structure (Fig. 2) from Eqn. 
Q). They correspond to additional hopping processes 
not taken into account in the simple model with only 
nearest-neighbor hopping. Comparing the Si- with the 
Ge-compound we may note a rough agreement with re- 
spect to the inter-chain hoppings. But, as can be also ex- 
pected from the different bonding angles 0, the in-chain 
hopping is clearly enhanced for BaCu2Gc207 in compar- 
ison to its Si-counterpart. 
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TABLE I: Collection of experimental data and theoretical 
model constants for BaCu2X207 (X=Ge,Si). 

We may use the derived TB parameters to construct a 
highly correlated one-band Hubbard model 

i,±z,a i,±x/y,o 

+ U ^2 n hl n i,l ( 2 ) 

i 

by identifying t z = t and using for t± the average value 
of t x and t y . This Hubbard- like Hamiltonian may be fur- 
ther reduced to a Heisenberg one H = (1/2) j JijSiSj, 
where the exchange couplings are connected with the TB 
parameters by Jy = AtfJU . Please note, that for the 
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purpose of estimating the inter-chain couplings, we ne- 
glected those terms in the Hamiltonian which are the 
most interesting in the present case, namely the DM term 
and the symmetric anisotropy. They are caused by the 
spin-orbit coupling ASL (A « 0.1 eV) and should be un- 
derstood to be added to the isotropic part. A pertur- 
bative estimate leads to an orbital moment oc X/sd (£d - 
crystal field splitting of d-levels) and a spin canting angle 
of 3° caused by the DM interaction oc A. The symmetric 
anisotropy oc A 2 compensates partly the DM term and it 
is just the Hund's rule coupling which leads to a "residual 
anisotropy" and to the two gaps in the magnetic excita- 
tion spectruniiii 

Experimental information about the nearest-neighbor 
exchange coupling J was derived from the maximum 
of the susceptibility curves^ using the Bonner-Fisher 
theory^ (see Table HJ. Using simply the same U for both 
compounds one would predict the ratio J Ge /J Sl = 1.84 
from the theoretical t = t z values, which is already in 
good agreement with the experimental ratio 1.9. One 
can also take the experimental J-values and the band- 
structure t's to derive an effective one-band Hubbard U 
of 4.8 cV (Gc) and 5.1 eV (Si) which seem to be quite 
reasonable (see Table HJ . In the given case the simple 
formula J = At 2 /U works much better than for Sr 2 Cu03 
or CuGeOs where considerable ferromagnetic corrections 
are necessary^ Due to the zigzag chains the hopping t 
parameters of Tableware much smaller than for Sr 2 Cu03 
which gives a better justification of the second order per- 
turbation theory. On the other hand, the Cu ions are 
more far away from each other than in CuGc03 which 
reduces the direct ferromagnetic exchange contributions. 

The theoretical inter-chain hopping and exchange pa- 
rameters are connected with quite important error bars 
(which is, however, not unusual for cuprate compounds). 
One possibility to derive the magnetic inter-chain cou- 
plings was discussed above and leads to J°® nd = 4t 2 L /£7 
(see Table . Another method uses the theory of Irkhin 
and Katanini^ which is an improved Schulz theory^ of 
weakly coupled spin chains. It was recently confirmed 
by detailed quantum Monte Carlo calculations^ In the 
theory of Irkhin and Katanin, the Neel temperature Tat 
is given by 



T N = kJj_z ± xoL(AJ/T N ) 

with the numerical constants k = 0.70, xo 
A = 5.8, the number of nearest neighbors z± 
the scaling function 
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2.1884, 
: 4, and 
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with B = 0.15. Applying and Q and using the mea- 
sured values of Tn^ xp and J exp we derived the empirical 
inter-chain exchange constants J^ mp collected in TableQ] 
Both methods lead to more or less identical inter-chain 
couplings for the two compounds. But the J^°- nd are 
roughly a factor of two larger than the empirical values 
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TABLE II: Distribution of magnetization in BaCu2X2C>7 with 
X=Ge,Si using FPLO within the LSDA and the LSDA+[/ 
schemes. Numerical values are in Bohr magneton. 



jemp _ s i m ii ar discrepancy had been alread y ob served 
for Sr 2 Cu0 3 , Ca 2 Cu0 3 and CuGe0 3 in Ref.lll where 
also possible origins of that discrepancy had been dis- 
cussed. Our present estimate gives only a rough idea 
about the absolute values of inter-chain couplings, but 
cannot answer whether these couplings would be ferro- 
or antiferromagnetic. Certainly, a more precise determi- 
nation of inter-chain couplings deserves more investiga- 
tions. 
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FIG. 3: (Color online) Total and partial density of states 
(DOS) of BaCu 2 Ge 2 07 from LSDA+C/ calculation. The thick 
black line is the total DOS whereas the thin green and blue 
lines are the DOS of all coppers and oxygens orbitals, respec- 
tively. 

By allowing an antiferromagnetic (AFM) arrange- 
ments of spins we found an insulating AFM solu- 
tion already within the local spin-density approximation 
(LSDA). However, the obtained gap values of about 0.5 
eV arc unrcalistically small in LSDA. Therefore, we per- 
formed LSDA+L7 calculations with the Slater-Coulomb 
parameters F 2 =8.6eV, F4=5.4eV (corresponding to J = 
1 eV) and two values for Fq = U = 3.7 eV and 5 eV. 
The former value produces gaps corresponding to the 
deep blue color of the samples (Fig. but the latter 
one is more close to parameters common for cuprates^i 
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FIG. 4: (Color online)Local magnetization distribution of 
BaCu2Ge2 07 in the be plane passing through the vertical 
chain axis c. Shown is the spin density in a logarithmic scale: 
yellow to red (blue) parts correspond to spin up [u] (down [d]) 
density, green parts are very low local density. Dashed lines 
represent the elementary cell. 



We used the same optimized basis orbitals like for the 
LSDA calculation and the "around mean field" version 
of LSDA+C/. The total energy differences between FM 
and AFM solution per magnetic ion AE lead to the near- 
est neighbor exchanges J = A.E/m2 of 46.6 (65.2) meV 
and 22.6 (45.5) meV for the Ge- and Si-compound with 
U = 5.0 (3.7) eV. 

The magnetization distribution on the different sites 
of the two compounds within LSDA and LSDA+J7 (for 
Fq = 3.7 eV having only small deviations to the Fq = 5 
eV solution) is shown in Tablc[np£ In the plaquctte, each 



Cu ion is surrounded by 4 oxygens: 01 and 02 are the 
side oxygens, and 03 the bridging oxygen between two 
plaquettes. As one can see, a non negligible part of the 
magnetization is situated at the oxygen sites, mainly at 
01 and 02. Compared to the LSDA results, the proper 
inclusion of Coulomb correlations leads to an increase of 
magnetization at the Cu sites since it stabilizes the AFM 
long range order by decreasing the quantum fluctuations. 

To see the magnetization distribution in more detail we 
calculated the spin density n-\ — in real space within 
the elementary cell. It is presented in Fig. 0|m a loga- 
rithmic scale. We neglect here the small deviations from 
a collinear spin arrangement, which is justified by the 
small spin canting angle of 1° . . .2°. One observes the 
characteristic form of 3d orbitals at the Cu sites and also 
the different oxygen sites. Quite remarkable are the dif- 
ferent magnetization directions in the two lobes of the 
bridging oxygen, which explains the small numbers for 
77io3 ■ A sign change of the magnetization at 03 is the- 
oretically possible by sp hybridization. On the contrary, 
the magnetization points into one direction at the side 
oxygens. 

Summarizing, we determined model parameters and 
predicted oxygen moments in BaCu 2 X 2 07 (X=Ge,Si). 
The absolute values of inter-chain couplings | Jj_ | are very 
similar for BaCu2Ge207 and BaCu2Si20y, but the near- 
est neighbor intra-chain J's are different. The larger J 
for the Ge-system makes it to be a more ideal model 
compound than BaCu2Si207. Our calculated oxygen mo- 
ments should be observable by 1 O nuclear magnetic res- 
onance, by //-spin rotation measurements, or by neutron 
scattering. They are similar (but slightly smaller) than 
those in Li2Cu02 which were recently detected^ 

We thank A. Stcpanov and H. Rosner for interesting 
discussions. 
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